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Abstract 

The turbulent transport of impurity or minority species, as for example Tritium, is investigated in drift- 
Alfven edge turbulence. The full effects of perpendicular and parallel convection are kept for the impurity 
species. The impurity density develops a granular structure with steep gradients and locally exceeds its 
initial values due to the compressibihty of the flow. An approximate decomposition of the impurity flux into 
a diffusive part and an effective convective part (characterized by a pinch velocity) is performed and a net 
inward pinch effect is recovered. The pinch velocity is explained in terms of Turbulent Equipartition fl] and 
is found to vary poloidally. The results show that impurity transport modeling needs to be two-dimensional, 
considering besides the radial direction also the strong poloidal variation in the transport coefficients. 

PACS numbers: 52.25.Gj, 52.35.Ra, 52.65.Kj 
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The transport properties of impurities are of great concern in magnetic fusion devices. In par- 
ticular impurities accumulate in the center of the plasma, where they are responsible for significant 
radiation losses. It is well established that turbulence is the key drive for plasma transport in the 
edge region and thus turbulence will also dominate the transport of impurities in that region. While 
in neutral fluids the transport of passive scalars is a very active field of research [2], it has in the 
context of plasma turbulence not yet found corresponding attention and measurements of impu- 
rity transport are most often interpreted within reduced ID transport models jsl 0]. In plasma 
turbulence the velocity field is in general compressible what makes the transport and mixing more 
complex than in incompressible flows [2]. This puts strong demands on the numerical modeling. 

Pinching of impurities has been observed as well as fast inwards transport of trace Tritium in 
JET Tritium puffs |5]. Especially in the edge region an inward convective flux (pinch) of impurites 
is found, which is proportional to the measured diffusion [3]. In the core where anomalous, turbu- 
lent transport is less important than in the edge the observed inward pinch of impurities seems to 
be in good agreement with neoclassical predictions based on trapped particle orbits. The Pfirsch- 
Schliiter impurity transport in the edge 1^, leads to an inward impurity pinch due to ion-impurity 
collisions. However, there is no explanation for the high inward pinch velocities needed to de- 
scribe impurity transport through in the turbulent edge region and it's scaling with the effective 
diffusion. 

Here we investigate the transport of impurities as a passive tracer field in electromagnetic edge 
turbulence described by a standard model of drift- Alfven turbulence flSyi- Flux tube geometry 
is used, with local slab-like coordinates {x,y,s) lldll . The following equations for the fluctuations 
in density n, potential (|) with associated vorticity CO = V^(|), current J and parallel ion velocity u 
arise in the usual drift-scaling: 

^ + {(^,a)} = ^(n)+V„7+^a)Vi(0, (la) 

^ + nsQ + «} = ^ (n - 4)) + V„ (7 - u) + A/nV^n, (lb) 

^ (pA„+A?7) +71 {(^ J} = V„ {nEQ + n-(^) -CJ, (Ic) 

e (^^ + {(^, "} j = -Vh {heq + n) . (Id) 

In addition to the cross-field advection, the parallel derivatives carry non-linearities entering 
through All, while the operator 3C represents effects of normal and geodesic curvature due to mag- 
netic field inhomogeneity with (Og = 2^, R being the tokamak major radius and L± the mean 



gradient length of the pressure profile, 
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The parallel current J is connected to the magnetic potential given by 7 = — V^A,,. The parameters 
reflect the competition between parallel and perpendicular dynamics, represented in the scale ratio 
£ = {qR/L±)^ where q is the safety factor. The electron parallel dynamics is controlled by the 
parameters 

p=— P=tt£^ C = 0.5\—-iJ = Vp, (2) 
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with Xg being the electron collision time and the factor 0.51 reflects the parallel resistivity Hllll . 
heq is an equilibrium density associated with corresponding neoclassical fields and currents. In 
detail the curvature operator on the equilibrium density balances the neoclassical Pfirsch-Schliiter 
current Jps 

^(«£e) + V„7ps = 0, (3) 
where the Pfirsch-Schliiter current is driven by the corresponding neoclassical potential 

V,(^P5)=C7p5- (4) 

To determine the evolution of the impurity density ni^p we assume that the contribution of the 
impurities to the gross plasma density o is negligible, i. e. = n,^o + ";mp ~ /^kO- In experiments 
the assumption n,.o ^ f^imp is possibly not always fulfilled, especially not during the initial phase of 
a Tritium gas puff, where locally in the edge region the Tritium migth contribute a significant part 
of overall plasma density. The bulk plasma is quasineutral with ng = hi q, allowing to regard either 
the ion- or electron continuity equation to determine the density evolution. We can, however, not 
use a corresponding relation for nimp- For cold impurities the drift velocity is given by the E x B- 
and the ion polarisation drift: 

M 2 

dtriimp = ^^-L ■ {nimpdt^ - nimpK {^) - V|| {riimpU) - fJimp"^ ±ni,np (5) 

Here we have introduced the relative mass of the impurities M = Mimp/Mi and Z indicates the 
charge state of the impurity ions. The fluctuating quantities (|) and u are input from the dynamical 
evolution of the turbulence and the total time derivative dt includes advection with the compress- 
ible E X B velocity. 

The dynamical equation for the impurity density thus differs from the dynamical equation for 
the density fluctuations Eq. (flbb . Typical simplifications originating in the distinction between a 



background and fluctuations are not possible to introduce for the impurity species. Finite inertia 
effects of the impurity ions enter through the ion-polarisation drift and finally all convection terms 
need to be preserved, while for example parallel convection is neglected in the dynamical equation 
for the density fluctuations. Note that the equation for the impurity species do not reflect the 
properties of the equilibrium as given by Eqs. @ and (El) as the impurities are not assumed to 
contribute significantly to the establishment of quasi-neutrality. 

Simulations were performed on a grid with 128 x 512 x 32 points and dimensions 64 x 256 x 2n 
in x,y,s. Standard parameters for the runs were /j = 5, q = 3, magnetic shear 5^= 1 (appearing 
only in the geometrical setup of the simulations [7] ), and qob = 0.05, with the viscosities set to 
l^a = l^n = 0.025, corresponding to typical edge parameters of large fusion devices. 

In the saturated turbulent state the equilibrium density gradient is weakly flattened. Strong 
deviations from the initially specified density gradient are, however, prevented by a feed-back 
mechanism using two damping layers at the inner and outer radial boundary. This forces the 
flux-surface averaged density to stay close to its initially specified equilibrium profile. 

To investigate impurity diffusion in fully developed, quasi-stationary turbulence we numeri- 
cally solve Eqs. (1) and let initial perturbations develop until a state of saturated, quasistationary 
turbulence is reached. The impurities are then released into the turbulence and convected around 
by the turbulence according to the evolving turbulent velocity field. The initial impurity density 
Hiinp is chosen as a radially localized Gaussian added to a constant impurity background density. 
For some runs the impurity density was additionally localized along the magnetic field lines, that 
is in coordinate s, to investigate the effect of parallel convective transport. Here we choose to 
investigate the behavior of massless impurities. The coupled system of bulk plasma turbulence 
and impurities is evolved until significant mixing of the impurities has been achieved and initial 
transient effects have decayed. For each parameter several runs are performed to increase the 
statistical significance of the results. 

From the transport of passive fields in compressible fluid turbulence it is well known that the 
passive quantity reveals a much flatter fluctuation spectrum than the turbulent energy spectrum, 
moreover the passive scalar tends to accumulate in the contracting regions of the turbulent flow 
field [2]. Thus, the impurity density after some time not only exhibits strong gradients, but locally 
the initial value of the impurity density can be exceeded. To be able to handle these effects within 
the given limited resolution the diffusive term in Eq. Q was chosen as ^i^p = 5/j„ 

A prominent feature of the impurity behavior is the weak parallel convective transport com- 
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pared to the radial turbulent transport. The reason is that the impurities are convected in the par- 
allel direction by the fluctuating parallel ion speed u which is small u ~ 0.01 compared to a radial 
velocity that is of order one. This is clearly observed in Fig. [H which shows the impurity density 
projected onto a geometrically poloidal cut. The projection roughly translates the variation in s to 
a poloidal variation, s = corresponds to the outboard midplane, 5 = 7r/2 to the upper, s = —Tl/2 
to the lower side and finally s = ±n corresponds to the high field side. The radial extend of the 
simulation domain has been stretched to allow better visualization. The initial impurity density 
is localized at the outboard midplane, corresponding to a parallel localization in the flux tube ge- 
ometry. No significant parallel flow of the impurity density is observed, while significant radial 
mixing occurs. Parallel compressional effects are however visible and arrange for finite passive 
density gradients at the high field side. Moreover an inward pinch effect is clearly observed at the 
outboard midplane. 

By starting from an initial impurity distribution that is homogeneous along s, this pinching 
velocity is seen to lead to a shift towards the torus axis of the impurity density compared to the 
initial condition (Figure |2l). 

For a more quantitative description of this behavior, the flux F of the impurity ion species can 
in lowest order be expressed by a diffusion coefficient D and a convective velocity V, which is 
associated to a pinch effect: 

Tyis) = -D{s)d^ < n >y W{s) <n>y . 

The turbulence is radially homogeneous and so there is no radial dependence of D and V . Aver- 
ages are taken along the periodic y direction. We obtain these values for each value of s along the 
magnetic field lines and thus at different poloidal positions. From a scatter plot of T{r) / < n>y 
versus < Inn >y the parameters D{s) and V{s) can be obtained. These are standard parameters 
used in modeling and in evaluation of transport experiments. The fitting procedure is carried out 
using a standard nonlinear least-squares Marquardt-Levenberg algorithm jol as implemented in 
the GNUPLOT software. Figure |3] shows such a typical scatter plot with a fitted linear relation 
between the two quantities indicating that while there are significant deviations from a linear re- 
lationship, the decomposition of the transport into D and V has some merit. The comparison of 
the evolution of the impurity density profile with the analytical evolution of the profile using the 
values for D and V, obtained from an analysis as shown in Figure |3] is depicted in Figure IH 

The poloidal dependence of diffusion and effective convection is rather strong and depicted in 
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Fig. 121 The effective convective velocity V{s) changes sign and is at the high field side directed 
outwards. This pinching velocity is due to normal curvature and can be consistently explained in 
the framework of Turbulent EquiPartition (TEP) III 13]. In the absence of parallel convection, 



finite mass effects and diffusion Eq. ^ has the following approximate Lagrangian invariant 

L{s) = Inriimp + (iiBCOs{s)x — GiBsin{s)y . (6) 

Its spatial homogenization on each drift plane < L{s) >y= const (s) by the turbulence indi- 
cates that at the outboard midplane (s = 0) the impurites are effectively convected radially inward 
leading to a radially gradient (< Inni^p >j°« const. — cOsJc) while at the high field side they are 
effectively convected radially outward (< Inntmp >y°^ const. + a)gjc). One should note that this 
effective inward or outward convection is not found as an average E x B velocity, but is miti- 
gated by the effect of homogenization of L under the action of the turbulence. The strength of the 
"pinch" effect is consequently proportional to the mixing properties of the turbulence and thus to 
lowest order has an additional variation that scales with the measured turbulent diffusion, so that 
V{s) —cos{s)D{s). The slight ballooning in the turbulence level thus causes the inward flow 
on the outboard midplane to be stronger than the effective outflow on the high-field side. Aver- 
aged over a flux surface and assuming poloidally constant impurity density a net impurity inflow 
results. The net pinch is directly proportional to the diffusion coefficient D in agreement with 
experimental observations Translated to dimensional values for typical large tokamak edge 
parameters we obtain D{s) oc 1 .5 — 2.0 m^/ s and V{s) +60 80 m/ s and an flux-surface aver- 
aged inward convection velocity of < V >= —0.4 m/s. Locally at the outboard midplane values 
ofV{s = 0)/D{s = 0) ~ — 40m^^ are found, in rough agreement with experimental values ^3]. 

Furthermore we observe a slight shift of the peak diffusion coefficient in the diamagnetic di- 

n 

rection. This effect is likely due to the spreading of the turbulence II15I1 in that direction leading 
to enhanced turbulence levels slightly upward from the outboard midplane and consequently to an 
up-down asymmetry. 

The strong peaking of the impurity density, the poloidal dependence of the transport coeffi- 
cients and the slow parallel diffusion of the impurities thus make it necessary to apply at least two 
dimensional modeling. Both a poloidally varying diffusion coefficient and effective convection 
velocity should be used in transport codes to describe impurity density evolution. The observed 
impurity pinch in the edge plasma region can be explained by turbulent equipartition without in- 
voking arguments from neoclassical transport theory. 
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FIG. 1: Impurity distribution projected onto a poloidal cross-section (note: radial dimension not to scale), 
white: high impurity density and dark: low impurity density. Left: initial distribution, localized along the 

magnetic field on the outboard midplane (at = . Right: after 25 time units corresponding to about lOO^us. 
Parallel transport is slow compared to radial transport. The inward pinch effect is clearly visible. 
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FIG. 2: Impurity distribution projected onto a poloidal cross-section (note: radial dimension not to scale). 
Left: initial distribution. Right: after 25 time units. The inward pinch effect on the outboard midplane and 
outward convective transport on the high field side (inboard midplane) is obvious. 
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FIG. 3: Typical scatter plot (at high field side) of the flux versus gradient with Unear fit. 
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Evolution of a Gaussian, Outboard Midplane 
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FIG. 4: Impurity density averaged along y at the outboard midplane and compared to the expected evolution 
of a Gaussian from the fitted coefficients D and V. 
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FIG. 5: Impurity diffusion D (top) and pinch velocity V over poloidal position angle with indicated error- 
bars. The observed asymmetry is due to turbulence spreading in the diamagnetic direction. 
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